Slab Load Controls Beneath the Alps on the Source-to-Sink Sedimentary Pathways in the Molasse Basin

The stratigraphic development of foreland basins has mainly been related to surface loading in the adjacent orogens, whereas the control of slab loads on these basins has received much less attention. This has also been the case for interpreting the relationships between the Oligocene to Micoene evolution of the European Alps and the North Alpine foreland basin or Molasse basin. In this trough, periods of rapid subsidence have generally been considered as a response to the growth of the Alpine topography, and thus to the construction of larger surface loads. However, such views conflict with observations where the surface growth in the Alps has been partly decoupled from the subsidence history in the basin. In addition, surface loads alone are not capable of explaining the contrasts in the stratigraphic development particularly between its central and eastern portions. Here, we present an alternative view on the evolution of the Molasse basin. We focus on the time interval between c. 30 and 15 Ma and relate the basin-scale development of this trough to the subduction processes, and thus to the development of slab loads beneath the European Alps. At 30 Ma, the western and central portions of this basin experienced a change from deep marine underfilled (Flysch stage) to overfilled terrestrial conditions (Molasse stage). During this time, however, a deep marine Flysch-type environment prevailed in the eastern part of the basin. This was also the final sedimentary sink as sediment was routed along the topographic axis from the western/central to the eastern part of this trough. We interpret the change from basin underfill to overfill in the western and central basin as a response to oceanic lithosphere slab-breakoff beneath the Central and Western Alps. This is considered to have resulted in a growth of the Alpine topography in these portions of the Alps, an increase in surface erosion and an augmentation in sediment supply to the basin, and thus in the observed change from basin underfill to overfill. In the eastern part of the basin, however, underfilled Flysch-type conditions prevailed until 20 Ma, and subsidence rates were higher than in the western and central parts. We interpret that high subsidence rates in the eastern Molasse occurred in response to slab loads beneath the Eastern Alps, where the subducted oceanic slab remained attached to the European plate and downwarped the plate in the East. Accordingly, in the central and western parts, the growth of the Alpine topography, the increase in sediment flux and the change from basin underfill to overfill most likely reflect the response to slab delamination beneath the Central Alps. In contrast, in the eastern part, the possibly subdued topography in the Eastern Alps, the low sediment flux and the maintenance of a deep marine Flysch-type basin records a situation where the oceanic slab was still attached to the European plate. The situation changed at 20 Ma, when the eastern part of the basin chronicled a change from deep marine (underfilled) to shallow marine and then terrestrial (overfilled conditions). During the same time, subsidence rates in the eastern basin decreased, deformation at the Alpine front came to a halt and sediment supply to the basin increased possibly in response to a growth of the topography in the Eastern Alps. This was also the time when the sediment routing in the basin axis changed from an east-directed sediment dispersal prior to 20 Ma, to a west-oriented sediment transport thereafter and thus to the opposite direction. We relate these changes to the occurrence of oceanic slab breakoff beneath the Eastern Alps, which most likely resulted in a rebound of the plate, a growth of the topography in the Eastern Alps and a larger sediment flux to the eastern portion of the basin. Beneath the Central and Western Alps, however, the continental lithosphere slab remained attached to the European plate, thereby resulting in a continued downwarping of the plate in its central and western portions. This plate downwarping beneath the central and western Molasse together with the rebound of the foreland plate in the East possibly explains the inversion of the drainage direction. We thus propose that slab loads beneath the Alps were presumably the most important drivers for the development of the Molasse basin at the basin scale

Local and regional minimum 1D models for earthquake location and data quality assessment in complex tectonic regions: application to Switzerland

One-dimensional (1D) velocity models are still widely used for computing earthquake locations at seismological centers or in regions where three-dimensional (3D) velocity models are not available due to the lack of data of sufficiently high quality. The concept of the minimum 1D model with appropriate station corrections provides a framework to compute initial hypocenter locations and seismic velocities for local earthquake tomography. Since a minimum 1D model represents a solution to the coupled hypocenter-velocity problem it also represents a suitable velocity model for earthquake location and data quality assessment, such as evaluating the consistency in assigning pre-defined weighting classes and average picking error. Nevertheless, the use of a simple 1D velocity structure in combination with station delays raises the question of how appropriate the minimum 1D model concept is when applied to complex tectonic regions with significant three-dimensional (3D) variations in seismic velocities. In this study we compute one regional minimum 1D model and three local minimum 1D models for selected subregions of the Swiss Alpine region, which exhibits a strongly varying Moho topography. We compare the regional and local minimum 1D models in terms of earthquake locations and data quality assessment to measure their performance. Our results show that the local minimum 1D models provide more realistic hypocenter locations and better data fits than a single model for the Alpine region. We attribute this to the fact that in a local minimum 1D model local and regional effects of the velocity structure can be better separated. Consequently, in tectonically complex regions, minimum 1D models should be computed in sub-regions defined by similar structure, if they are used for earthquake location and data quality assessmen

Combining controlled-source seismology and local earthquake tomography to derive a 3-D crustal model of the western Alpine region

We present a newly developed approach of combining controlled-source seismology (CSS) and local earthquake tomography (LET) data to obtain a new 3-D crustal model of the western Alpine region. Our approach combines either data by taking into account the strengths of the individual seismic methods. Our western Alpine 3-D model is primarily based on a well-defined Moho, constrained by CSS and LET data, and includes smooth lateral variations in seismic velocities mainly constrained by LET data, but locally also by CSS data. The consistent combination of results from the two different seismic methods is feasible due to LET Moho elements, as defined by characteristic P-wave velocities and their uncertainty estimates. These uncertainty estimates are based on values of the diagonal element of the resolution matrix, absolute P-wave velocities that are typical for crust and mantle and a specific velocity gradient across the Moho discontinuity. Finally, our definition of LET Moho elements and their uncertainties is validated by comparisons of highest quality Moho results from both methods coinciding in 353 localities. Our model clearly shows three Moho surfaces, being Europe, Adria and Liguria as well as major tectonic structures like suture zones and the high-velocity Ivrea body. In general, it is in a good agreement with previous studies. The biggest differences occur along plate boundaries, where the strong lateral velocity variations are best resolved by LET. Due to the larger number of available Moho reflector elements a more accurate definition of plate boundaries at Moho level is possible and, therefore, new insights in deep lithosphere structures of the Alpine collision zone can be expected. Furthermore, our new 3-D crustal model directly includes a 3-D migrated image of the Ivrea bod

High-resolution body wave tomography beneath the SVEKALAPKO array — II. Anomalous upper mantle structure beneath the central Baltic Shield

A number of different geodynamic models have been proposed to explain the early tectonic evolution of the Baltic Shield. To provide additional geophysical constraints on these models, we performed a teleseismic tomography traveltime inversion for the central part of the Baltic Shield. The SVEKALAPKO project is focused on the investigation of the lithosphere—asthenosphere structure down to 400 km depth under central Fennoscandia (Baltic Shield). A total of 143 stations were deployed including 15 permanent stations from the Finnish seismic network. The temporal network was composed of 40 broad-band and 88 short-period instruments distributed in a rectangular array of 1000 km by 900 km from 1998 August to 1999 May. The results are based on a non-linear teleseismic tomography algorithm. They reveal significant P-velocity variations (up to 4 per cent) throughout the SVEKALAPKO array. The most prominent feature is a positive anomaly that can be followed down to 250 km depth beneath the centre of the array. We interpret this anomaly as the signature of the tectosphere (Jordan 1978) beneath the Fennoscandian Shield. It correlates spatially with an anomalous high-velocity lower crust. Other shallow (crustal) anomalies can be correlated with magmatic events surrounding this nucleus of high velocity. Comparison of images before and after correction by crustal structure proves that this methodology yields solid and coherent tomographic results. Further observations of relative P traveltime residuals from six teleseismic events with different azimuths show delay variations of ±2.0 s between stations located in the North German basin and stations on the Svecofennian Shiel

High-resolution body wave tomography beneath the SVEKALAPKO array: I. A priori three-dimensional crustal model and associated traveltime effects on teleseismic wave fronts

Assessment of contributions from shallow lithosphere to teleseismic wave front distortion is a prerequisite for high-resolution regional teleseismic tomography. Several methods have been proposed in the past for the correction of these effects, e.g. application of station correction terms. We propose an approach that is independent of the subsequent inversion and uses the available a priori knowledge of the crustal structure to calculate crustal traveltime effects of teleseismic wave fronts. Our approach involves the construction of a 3-D crustal model based on controlled source seismology data and calculation of the associated traveltime anomalies for incoming teleseismic wave fronts. The model for central Fennoscandia shows a maximum crustal thickness of 64 km and includes a high-velocity lower crust as derived for parts of the study area by previous authors. Traveltimes calculated using finite differences for teleseismic waves travelling through this crustal model are compared with those from the standard reference model IASP91 and the residuals are used to correct observed teleseismic arrival times at the SVEKALAPKO array. To test the performance of this approach, in a second part of the study a synthetic traveltime data set is obtained by tracing wave fronts through a mantle structure with known velocity anomalies and the 3-D crustal model. This data set is inverted with and without correction for crustal effects. The 3-D crustal effects alone with a homogeneous mantle are also inverted and the results show that the crustal effects propagate down to 450 km. The comparison of the inversion results demonstrates the need to apply appropriate 3-D crustal corrections in high-resolution regional tomography for upper-mantle structure beneath the Baltic Shiel

Induced seismicity during the construction of the Gotthard Base Tunnel, Switzerland: hypocenter locations and source dimensions

A series of 112 earthquakes was recorded between October 2005 and August 2007 during the excavation of the MFS Faido, the southernmost access point of the new Gotthard Base Tunnel. Earthquakes were recorded at a dense network of 11 stations, including 2 stations in the tunnel. Local magnitudes computed from Wood-Anderson-filtered horizontal component seismograms ranged from −1.0 to 2.4; the largest earthquake was strongly felt at the surface and caused considerable damage in the tunnel. Hypocenter locations obtained routinely using a regional 3-D P-wave velocity model and a constant Vp/Vs ratio 1.71 were about 2km below the tunnel. The use of seismic velocities calibrated from a shot in the tunnel revealed that routinely obtained hypocenter locations were systematically biased to greater depth and are now relocated to be on the tunnel level. Relocation of the shot using these calibrated velocities yields a location accuracy of 25m in longitude, 70m in latitude, and 250m in focal depth. Double-difference relative relocations of two clusters with highly similar waveforms showed a NW-SE striking trend that is consistent with the strike of mapped faults in the MFS Faido. Source dimensions computed using the quasidynamic model of Madariaga (Bull Seismo Soc Am 66(3):639-666, 1976) range from 50 to 170m. Overlapping source dimensions for earthquakes within the two main clusters suggests that the same fault patch was ruptured repeatedly. The observed seismicity was likely caused by stress redistribution due to the excavation work in the MFS Faid

High-resolution teleseismic tomography of upper-mantle structure using an a priori three-dimensional crustal model

The effect of an a priori known 3-D crustal model in teleseismic tomography of upper-mantle structure is investigated. We developed a 3-D crustal P-wave velocity model for the greater Alpine region, encompassing the central and western Alps and the northern Apennines, to estimate the crustal contribution to teleseismic traveltimes. The model is constructed by comparative use of published information from active and passive seismic surveys. The model components are chosen to represent the present large-scale Alpine crustal structure and for their significant effect on the propagation of seismic wavefields. They are first-order structures such as the crust-mantle boundary, sedimentary basins and the high-velocity Ivrea body. Teleseismic traveltime residuals are calculated for a realistic distribution of azimuths and distances by coupling a finite-difference technique to the IASP91 traveltime tables. Residuals are produced for a synthetic upper-mantle model featuring two slab structures and the 3-D crustal model on top of it. The crustal model produces traveltime residuals in the range between −0.7 and 1.5 s that vary strongly as a function of backazimuth and epicentral distance. We find that the non-linear inversion of the synthetic residuals without correcting for the 3-D crustal structure erroneously maps the crustal anomalies into the upper mantle. Correction of the residuals for crustal structure before inversion properly recovers the synthetic slab structures placed in the upper mantle. We conclude that with the increasing amount of high-quality seismic traveltime data, correction for near-surface structure is essential for increasing resolution in tomographic images of upper-mantle structur

Seismic moulin tremor

Through glacial moulins, meltwater is routed from the glacier surface to its base. Moulins are a main feature feeding subglacial drainage systems and thus influencing basal motion and ice dynamics, but their geometry remains poorly known. Here we show that analysis of the seismic wavefield generated by water falling into a moulin can help constrain its geometry. We present modeling results of hour-long seimic tremors emitted from a vertical moulin shaft, observed with a seismometer array installed at the surface of the Greenland Ice Sheet. The tremor was triggered when the moulin water level exceeded a certain height, which we associate with the threshold for the waterfall to hit directly the surface of the moulin water column. The amplitude of the tremor signal changed over each tremor episode, in close relation to the amount of inflowing water. The tremor spectrum features multiple prominent peaks, whose characteristic frequencies are distributed like the resonant modes of a semiopen organ pipe and were found to depend on the moulin water level, consistent with a source composed of resonant tube waves (water pressure waves coupled to elastic deformation of the moulin walls) along the water-filled moulin pipe. Analysis of surface particle motions lends further support to this interpretation. The seismic wavefield was modeled as a superposition of sustained wave radiation by pressure sources on the side walls and at the bottom of the moulin. The former was found to dominate the wave field at close distance and the latter at large distance to the moulin

High-precision earthquake locations in Switzerland using regional secondary arrivals in a 3-D velocity model

We present a new approach to relocate earthquakes in the greater western Alpine region using main crustal phases (Pg, Pn, PmP) that takes advantage of recent developments in P-wave velocity models and modelling of the Moho topography in the region, as well as the ability to track reflected and refracted phases in three-dimensional (3-D) heterogeneous media. Our approach includes a new 3-D P-wave velocity model for Switzerland and surrounding regions that combines a first-order Moho discontinuity based on local earthquake tomography (LET) and controlled-source seismology (CSS) information and 3-D seismic velocity information based on LET. Traveltimes for the main crustal phases (Pg, Pn, PmP) are computed using a fast marching method. We use a non-linear, probabilistic approach to relocate earthquakes that has been extended to include the use of secondary phases. We validate our approach using synthetic data, which was computed for a real earthquake and different combinations of available phases (Pg, Pn, PmP). We also applied our approach to relocate four selected earthquakes, two shallow and two deep crustal events in the northern Alpine foreland, for which independent information (ground truth information) on their focal depths exist. Our results demonstrate that the precision and accuracy of focal depth estimates can be greatly improved if secondary phases are used. This gain is a combined effect of an improved range of take-off angles and the use of differential traveltimes between first and secondary arriving phases. Our results also show that reliable information on the Moho depth is crucial to obtain accurate focal depths, if Pn or PmP phases are used in the relocation process. Finally, our approach demonstrates that proper identification of the main crustal phases in combination with an appropriate model parametrization in the forward solver will significantly improve earthquake location

3D crustal structure of the Eastern Alpine region from ambient noise tomography

The tectonic evolution of the European Eastern Alps within the Alpine orogeny is still under debate. Open questions include: the link between surface, crustal and mantle structures; the nature of the Moho gap between the two plates; the relationship between the Alps, the adjacent foreland basin and the Bohemian Massif lithospheric blocks. We collected one year of continuous data recorded by ~250 broadband seismic stations –55 of which installed within the EASI AlpArray complementary experiment– in the Eastern Alpine region. Exploiting surface wave group velocity from seismic ambient noise, we obtained an high-resolution 3D S-wave crustal model of the area.The Rayleigh-wave group-velocity from 3 s to 35 s are inverted to obtain 2-D group velocity maps with a resolution of ~15 km. From these maps, we determine a set of 1D velocity models via a Neighborhood Algorithm, resulting in a new 3D model of S-wave velocity with associated uncertainties. The vertical parameterization is a 3-layer crust with the velocity properties in each layer described by a gradient. Our final model finds high correlation with specific geological features in the Eastern Alps up to 20 km depth, the deep structure of the Molasse basin and important variations of crustal thickness and velocities as a result of the Alpine orogeny post-collisional evolution. The strength of our new information relies on the absolute S-wave crustal velocity and the velocity gradient unambiguously sampled along the Moho, only limited by the amount and quality distribution of the data available